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Abstract. The Properties of the Mott insulator to superfluid phase transition are obtained through the 
fcrmionic approximation in the Jaynes-Cummings-Hubbard model on linear, square, SC, FCC, and BCC 
Bravais lattices, For varying excitation number and atom-cavity frequency detuning. We find that the Mott 
lobes and the critical hopping are not scalable only for the FCC lattice. At the large excitation number 
regime, the critical hopping is scalable for all the lattices and it does not depend on the detuning. 
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Approximation. When the hopping term is non-zero the 
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1 Cold-atom systems have been regarded as efficient Simula- , ,. , , . . n r^, -..a, n , 

■ J ° solution becomes non trivial, the difficulty to rind ana- 
tors of quantum many-body physics 1 1 2 due to its ease of , ,. , , ,. c , , , c j i , 

n J J !— ^ ! lytical solutions lor the model has forced researchers to 

controllability. Research involving ultracold bosonic sys- . . , . , , r i i- 

resort to approximation or numerical methods lor dealing 
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^ — 1 . terns has brought lot of interest in the subject 13141 . Exper- ... , r™ n ,, . , , rm u 

° J r with such problems [13]. Recently, Mermg et ai. [14] have 

' ' iments with optical lattices in three dimensions have re- , , 

• • ■ proposed an approach m which spin operators are mapped 

L 1 "^ 1 vealed the superfluid-Mott insulator (MI-SF) phase tran- „ , „ . ,, , . „ 

■ v ' to the fcrmionic ones, hence allowing the application of a 

h3 ■ sition IS] . Such phenomena have been studied in the frame- _ . „ , , . , TT 
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workofthe Jaynes-Cummings-Hubbard Model (JCHM) EEEE]. , , , . , . , , , , 

J ° v / . uL-pcndcnt ones, which arc associated to each momentum 

The J CHM has been a widely used tool for the investi- yalue The great advantage of this method is the simp iicity 

gation of many-body systems describing the interplay be- in which physical quantit i es are found; such as the ener- 

tween the atom-cavity coupling and the inter-cavity hop- gigs &nd ^ chemical potential. The approach presented 

ping of photons QUI- In the absence of the hopping term, by Mering ^ al includes all classes of Bravais structures . 
the model reduces to the Jaynes-Cummings Model |llll2j 
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However, they only derived the results for one-dimensional When t = 0, the Hamiltonian (JXJ) is decoupled into 
lattices. L independent Jaynes-Cummings model Hamiltonians. In 
Despite the experiments involving optical structures in this case the system has well-known eigenstates |16) . For 
three dimensions, only a few theoretical results regarding i ^ 0, the atoms becomes coupled thus increasing the 
these lattices in dimensions greater than one is available in complexity of the solution, Since we can not write the 
the current literature |15j . It lacks a systematic presenta- eigenstates of the whole system as a direct product of 
tion of the phase diagram of typical Bravais lattices. This single-cavity eigenstates. As discussed in the introduction, 
is the purpose of the present paper. Here, we investigate an appropriate approach is the fermionic approximation 
the JCHM for different Bravais lattices in one, two and recently introduced by Mering et al. [14] . 
three dimensions and analyse the influence of the topol- 
ogy on the MI-SF phase transition. We use the approach 3 The Fermionic Approximation 
introduced by Mering et al. [14] . 



The paper is organized as follows. In Section II we in- 
troduce the JCHM. In Section III we present the fermionic 
approximation. Results are presented in Section IV. Fi- 
nally, IN Section V, we summarize our main results and 



The fermionic approximation consists in replacing the spin 
operators by fermionic ones, i.e., a + (o" — ) IS replaced by 
(c). In this framework, we can rewrite Hamiltonian (TXJ) 
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This approximation allows to solve the model exactly, by 



rp, T<"itji\/r u ma • c i j-j-- t t 4- • • means of a Fourier transformation. For t — 0, spin and 
I he JCHM hamiltonian for a lattice of L atoms is given ' * 

k _ ^ fermionic operators are equivalent and then the approxi- 
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mation becomes exact. Therefore, for small values of t this 
approach turns out to be very accurate when dealing with 
the JCHM [H]. 



where a ± = a x ± ia v and a x , y , z are the usual Pauli matri- Now we a ^ a Fourier transform to the fermionic and 
ces, aj (at) is the annihilation (creation) operator of the bosomc operators as 

1 k Ft ■ 1 k Ft ■ 

light mode at the jth atom, uj is the light mode frequency, = _ ^ e _27ri_ c i a fc , Cj = —j= ^ e~ 27ri_ ^c fc , (3) 

and e the atomic transition frequency. The light-atom cou- 



then the Hamiltonian can be written as 

pling is represented by g, t is the hopping integral, and (ij) 

H = 

denotes pairs of nearest-neighbour atoms on the lattice. 



^ [wfcaj^afe + g{a\c k + a k c\) + ec\c k , (4) 
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where u k = <->J-v k and v k is the dispersion relation of the 
Bravais lattice. 

The Hamiltonian (0| corresponds to a sum of L inde- 
pendent Hamiltonians H k {H = ^2 k H k ), where each of 
them is associated with a particular momentum k. The 
ground-state energy of H k is given by 
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(5) 



where the superscript denotes the excitation number and 
A = e — uj is the detuning between atomic trantition and 
light frequencies. Notice that h k commutes with H k . For a 
total excitation number N (N = ^2 k n k ) we have a partic- 
ular configuration {n kl ,n k2 ,n k3 ^ ...} that minimizes the 
total ground-state energy, J2 k ^ k h - For t <IC 1, it is easy 
to see that this configuration is {n, n, n, ...} corresponding 
to the Mott insulator state. By increasing t, a quantum 
phase transition takes phace and the system is driven to 
a superfluid state. Since n is constant, the phase bound- 
aries of the Mott lobes are n dependent. Thus, the nth 
Mott lobe is obtained through an analysis of the particles 
chemical potential, [i + = E?, — E^, and the hole one, 
fi~ = E^ — E^ -1 |14I17] . where k' and k are, respectively, 
the values that minimize and maximize these potentials. 
For /i + = the Mott lobe is closed at the critical hop- 
ping, t c , hence describing the MI-SF transition. 

4 Results 

In order to analyse the influence of Bravais lattices topol- 
ogy on the MI-SF transition, we study the one-dimensional 
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Fig. 1. First three Mott lobes for different lattices. Inside 
the lobes we have a Mott insulator state, while outside the 
system is in a superfluid state. 

(ID), square (SQ), simple cubic (SC), body-centered cu- 
bic (BCC) and face-centered cubic (FCC) lattices. The 
dispersion relations v k are, respectively, given by (18) 



,(i-D) _ 



-2t cos(fca) 



v ^k^k ~ — 2t[cos(fc x a) + cos(fc i/ a)] 
^kfky k z = — 2t[cos(fc K a) + cos(k y a) + cos(k z a)] 
vj?^ C \ z = — 8t[cos(k x a) cos(k y a) cos(fc z a)] 



(6) 

(7) 
(8) 
(9) 



and 



v\ k ! — — 4t[cos(k x a) cos(k y a) + cos(k x a) cos(k z a) + 



cos(k y a) cos(fc 2 a)], 



(10) 



where a is the lattice constant. For each structure we found 
the momentum vectors that maximizes the hole chemical 
potentials and minimizes the particle ones in order to ob- 
tain u~ , u + , and consequently the Mott phase boundary. 

Figure [T] shows the first three Mott lobes for uj/g — 1 
and A — 0, for the considered lattices. We see that as 
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the number of lattice neighbors increases, the MI phase 
region decreases. This result is expected since the proba- 
bility of ohoton tunneling increases for higher number of 
nearest-neighbors. We observe that for <i-dimensional hy- 
percubic lattices the lobes can be rescaled by = t/d 
causing the collapse into a single curve. The shape of the 
Mott lobes are equal for bipartite lattices, where a bipar- 
tite structure is such one that we can decompose into two 
substructures, with all nearest-neighbour sites shared be- 
tween each other. The FCC lattice is non-bipartite, Hence 
displaying a different behaviour. Thus we propose to make 
a detailed comparative analysis between the FCC and SC 
lattices representing, respectively, the non-bipartite and 
bipartite classes. 

The first Mott lobe (n = 1) on the SC and FCC lat- 
tice is show in figure [5] for typical detuning values. We see 
that both lattices have a similar MFSF phase transition 
frame. However, the FCC always has a smaller MI phase 
region. For both lattices, the MI phase region decreases for 
increasing detuning. Figure [3] shows the critical hopping 
t c as a function of the detuning. While at the first lobe, 
t c decreases when A increases, in the other lobes (n > 1) 
we observe that the critical hopping reaches a maximum 
at A — A rn . For this particular detuning value, when n 
increases the critical hopping decreases, as predicted in 
[T3] . This behavior is present in both structures. Figure S] 
shows the critical hopping as a function of n for A = 0. 
The properties expressed by both lattices are again qual- 
itatively equivalent where there is only one gap between 
the two curves. Figure [5] presents the detuning values cor- 




(a) 



-0.2 



-0.4 
-0.6 — 



-0.8 



FCC 



A/g=-1 

A/g=0 

A/g=1 

A/g=2 



3 -to 

A 

w -1.2 
-1.4 
-1.6 
-1.8 



0.00 



0.02 



0.04 0.06 

t/g 



00 



0.08 



0.10 



Fig. 2. Mott lobe for n = 1 and typical values of detuning 
for (a) SC and (b) FCC lattices. 

responding to the maximal critical hopping as a function 
of n. We can observe that, except for the n = 1 case, 
where the t c maximum is associated with the detuning 
minimum, The critical hopping correspondent detuning is 
approximately null. 

By performing an asymptotic approach to large exci- 
tation number, n» 1, we can find the dominant term of 
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Fig. 3. Relationship between critical hopping and detun- 
ing on the (a) SC and (b) FCC lattices for varying exci- 
tation number. 

the critical hopping which is given by 
9 
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Fig. 4. Critical hopping for the SC and FCC lattices in 
terms of n for null detuning. 




t r . = 



+ 0{n~ 6 ) 



(11) 



where we obtain d — 4 for the FCC and BCC lattices 
where it corresponds to the hypercubic lattices dimen- 
tions, i. e., 1 for linear, 2 for square, and 3 for the SC 
lattice. It is important to emphasize that at the large-n 
regime, the detuning and the topology class (bipartite or 



Fig. 5. Detuning versus n for maximal critical hopping 
on the SC and FCC lattices. 

non-bipartite) influence on the critical hopping t c is sup- 
pressed. Figure [5] confirms the prediction of the equation 
(JTTJ) . It shows that the exact results are in excellent agree- 
ment with the asymptotic ones for n > 4. 

5 Conclusions 

We have studied the properties of the MI-SF phase tran- 
sition for the Jaynes-Cummings-Hubbard model over sev- 
eral Bravais lattices by means of the fermionic approxima- 
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in n which does not depend on A and topology associated 
parameter can be rescaled for different classes through a 
multiplication of t c by an effective parameter that cor- 
responds to the dimension, for hypercubic, and to 4, for 
BCC and FCC lattices. 
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Fig. 6. Critical hopping versus n for various lattices at 
the large n regime. The symbols are related to numeri- 
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